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Abstract. The problem of nonlinear filtering of a random field observed in the presence of a noise, 
modeled by a persistent fractional Brownian sheet of Hurst index {Hi, H2) with 0.5 < Hi,H2 < 1, 
is studied and a suitable version of the Bayes' formula for the optimal filter is obtained. Two types 
of spatial "fractional" analogues of the Duncan-Mortensen-Zakai equation are also derived: one 
tracks evolution of the unnormalizcd optimal filter along an arbitrary "monotone increasing" (in 
the sense of partial ordering in M^) one-dimensional curve in the plane, while the other describes 
dynamics of the filter along paths that are truly two-dimensional. Although the paper deals with 
the two-dimensional parameter space, the presented approach and results extend to d-parameter 
random fields with arbitrary d > 3. 
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in 

^ ! 1 Introduction 

. An important estimation problem, arising in many engineering and physical systems evolv- 

I ing in time and space, is that of recovering a signal {Xt, t E T) from an observed noisy 

' nonlinear functional of the signal, represented by a process (It, t G T). In the classical 

mathematical filtering framework, one has T = [0, 00) or T = [0, T], with t interpreted as 
^ I "time", and the problem then is to characterize the conditional distribution of Xt given 

\—{ ' the observation u-field J-^ = a{Ys,0 < s < t}, where the latter represents information 

supplied by the noisy observation process from time up to time t. However, there is 
a number of interesting applications, arising, for example, in connection with denoising 
of images and video-streams, where the parameter space T has to be multidimensional, 
which renders the classical theory of nonlinear filtering inapplicable. The latter observa- 
tion stems directly from the fact that, unlike M which permits perfect ordering, there is 
only partial ordering available in with d > 2, thus, on the one hand, use of the mul- 
tiparameter martingale theory in the underlying analysis is required, while, on the other 
hand, evolution of the optimal filter for d-parameter random fields can be studied along 
arbitrary ^-dimensional "monotone increasing" paths with 1 < i < d. 

To extend the classical one-parameter nonlinear filtering theory to the multiparameter 
spatial filtering case, it is natural to start with the following observation model for a 
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random field {Xt : t e T), with T = [0, Ti] x • • • x [0,Td], corrupted by an additive 
multiparameter observation noise M = {J\ft :teT): 

= / •••/ h{X^,^^,,,^,^))dsi...dsd+M(t,,...,t^), {ti,...,td) €T, (1) 
Jo Jo 

where his a (suitably integrable) nonlinear function of the "signal" of interest X. Consider 
the observation cr-field 

^l = ^ltr,...,t,)--=(^{Ys:0^s^t}, (2) 

with = (0, . . . , 0) G T and where, for all t = (ti, . . . , to) and s = (si, . . . , s^) in T, we 

put s ~< t whenever Sj < ti for all i = 1, . . . ,d. Then the aim of the filtering theory is 
to describe the conditional distribution of the true "signal" of interest X at "location" 
t, given the observation sigma-field J^Y o'^; equivalently, one can study the dynamics of 
E{F{Xt) I J-'Y) for a sufficiently rich class of test functions F. 

Interestingly, even in the case when observation noise in (1) is a standard two- 
parameter Wiener sheet, assumed to be independent of X (d = 2 here), h is square- 
integrable function and the signal is known to have a semimartingale structure, the actual 
derivation of evolution equations satisfied by the optimal filter is somewhat non-trivial, 
owing to the fact that the multiparameter martingale theory is significantly more com- 
plicated than the classical one and many "standard" martingale tools available in the 
one-parameter case are no longer applicable in the multiparameter setting. This formula- 
tion of the nonlinear spatial filtering problem (i.e. with standard Wiener sheet observation 
noise) has been studied in [1] and [2], where several types of stochastic partial differential 
equations governing the unnormalized optimal filter were obtained. However, for the case 
of other types of continuous multiparameter random fields J\f driving the observation field 
Y, no mathematical theory of optimal nonlinear filtering currently exists. 

The goal of the present paper is to study the above problem of nonlinear filtering 
of semimartingale random fields (with d > 2) but in the presence of a long-memory 
observation noise M, where the latter is modelled by a persistent fractional Brownian 
sheet. 

The paper is organized as follows. The remainder of Section 1 is devoted to two topics 
of interest: i) preliminaries on multiparameter martingales, which will be useful to us in 
Section 2; ii) properties of fractional Brownian sheet, plus a number of relevant results 
from fractional calculus. Section 2 presents the main theorems of the paper. Namely, 
as a first step, an appropriate spatial version of the "fractional" Bayes' formula is ob- 
tained. Next, a stochastic evolution equation for the unnormalized optimal filter along 
a one-dimensional monotone "increasing" path is derived. Finally, a stochastic evolution 
equation describing dynamics of the optimal filter along proper two-dimensional paths in 
the plane is also presented. Unlike what happens in the case of standard Wiener sheet 
observation noise, the latter two evolution equations cannot, strictly speaking, be inter- 
preted as measure- valued stochastic partial differential equations due to the effects of long 
memory, but they certainly represent the "fractional" multiparameter analogues of the 
classical Duncan-Mortensen-Zakai filtering equations. Lastly some concluding remarks 
are given in Section 3. 
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For the sake of brevity and notational convenience, from now on we will restrict our 
attention to the case of two-dimensional parameter space T, since, although analogous 
techniques and results can certainly be developed for the higher-dimensional cases, the 
latter lead to a larger number of terms in evolution equations and more cumbersome 
notation throughout the derivations. 

1.1 Two-pcirameter martingale theory 

First let us define the usual partial ordering -< in the positive quadrant M^, along with 
the following relations and operations: for arbitrary a = (ai,a2) and b = (61,62) in M^, 
a ^ 6 if and only if ai < bi and 02 < 62; a b if and only if ai < 61 and 02 < 62; 
a X 6 if and only if oi < 61 and 02 > 62; a A 6 := (min(ai, 61), min(o2, 62)); 
a V 6 := (max(ai, 61), max(a2, 62)); a 6 := (ai, 62). 

Given a random field X with a parameter set M^, define its "increment" over an 
arbitrary rectangle {z,z'] := {x = {xi,X2) G : xi € (2:1,2:^], X2 & (-2^2,-^2]}' where 
z -< z', by 

Next, for a complete probability space {^,!F,P), let {TziZ G M+} be a family of sub-cr- 
ficlds of T satisfying the following properties: 

(Fl) \iz<z' then C J^^/; 

(F2) contains all P-null sets of T] 

(F3) for each z G R^, Tz = ^z<<z' ^z'\ 

{FA) for each z = (^1,22) £ R+j and are conditionally independent given J^zi where 
J^l and J^l are defined by 

■■= \l ^{z^,t) = ^{ U -^(^i-t)} ■■= V ^(^.^2) = ^{ U -^(^.^a)}- 

t>0 t6]R+ s>0 selR+ 

Note that condition (F4) is equivalent to the following condition {FA'): for all bounded 
random variables X and all z G M^, 

E{X1.F4 = E{E{X\J^l}\J^^} =E{E{X\J^^}\J^l} a.s. 

Moreover, for an arbitrary random field X with independent increments, i.e. such that 
X{Ai), . . . , X{An) are independent for all disjoint rectangles Ai, . . . , An C M^, the natural 
filtration generated by X over rectangles satisfies condition (F4) (see [3]), i.e. JF^ := 
(7{X(^) : A ^ z} has property (F4), where we say that ^ -< z if x ^ 2; for all a; G A. 

Definition 1.1 Let {J^z)zeR'^ filtration satisfying (F1)-(F4). The process X = 

{Xz,z G M^} is called a two-parameter martingale with respect to {Tz) if: i) for each 
z G M^, Xz is adapted to J^z o,nd integrable; and ii) for each z -< z' , E(^Xzi\J^z) = ^z o-s- 
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Definition 1.2 Let X = {X^ : z G R^} be a process such that X^ is integrable for all 
z G and let filtration {J-'z)z^m^ satisfy (F1)-(F4). Then 

(a) X is called a weak martingale with respect to [J-z) if: 

(i) X is adapted to {Tz), o.nd (ii) E{X((z, z'])!^^} = a.s. for all z z' . 

(b) X is called an i-martingale (i = 1,2) with respect to (J^z) if' 

(i) Xz is J=i-adapted, and (ii) E{X{{z, z'])\J^^} = a.s. for all z -<-< z' . 

(c) X is called a strong martingale with respect to {J^z) if: 

(i) X is adapted to {J^z), 

(ii) X vanishes on the axes (i.e. ^(0,^2) ~ ^ and X(^z-ifi) = a.s. for all zi,Z2 & ^+), 
and (Hi) E{X{{z, z'])\Tj\/ T^} = a.s. for all z z' . 

Note that a martingale is both a 1- and a 2-martingale. The converse wih also hold (i.e. 
if X is both a 1- and a 2-martingale, then X is a two-parameter martingale), provided 
that o))-^(2^ 0)'-^! ^ ^+J' {"^^(0,22)' -^(022)' ^2 £ R+l are both martingales. Also, 

clearly, any martingale is a weak martingale and any strong martingale is a martingale. 

Let us say that a process {-X'^} is right- continuous if for a.e. lo, lim2,'_>2; Xzi{lo) = Xz{uj) 

z-<z' 

for all z G M^, and that it has left limits if, for a.e. \mi^i_^^ Xzi{oj) exists for all 

z'-«z 

zG(R+\{0})2. 

Definition 1.3 Given filtration {J^z) satisfying properties {Fl) — (-F4), a process X = 
{Xz,z G R^} is called an increasing process if: (i) X is right- continuous and adapted to 
(J^z), (ii) ^z = O'-s. on the axes, and (Hi) X{A) > for every rectangle A C M^. 

For our purposes, it will be sufficient to work with a bounded subset T = [0,ri] x [0, 

of instead of all of R^. Let us fix an arbitrary T = (Ti,T2) G R^, and, for p > 1, 
define MP(Y) to be the class of all right-continuous martingales M = {Mz,z ~< T} such 
that Mz = a.s. on the axes and E\Mz\p < 00 for ah z € T. Let M^{T) and M^g{T) 
denote respectively the class of continuous and the class of strong martingales in MP{T). 

The following result highlights some of the fundamental differences between the clas- 
sical one-parameter martingale theory and its multiparameter analogue. As shown by 
Cairoli and Walsh in [3] in the two-parameter case, for an arbitrary martingale M G 
M^{T), there exists an increasing process A = {Az,z G T} such that {M| — Az,z G T} 
is a weak martingale. However, such an increasing process A need not be unique even in 
the case of a strong martingale M. Nor can one in general guarantee the existence of an 
increasing process A such that {M^ — 2: G T} is a regular two-parameter martingale. 
Thus, we will agree to denote by (M) = {{M)z,z G T} any increasing process A such that 

— ^ is a weak martingale. Some refinements of the above weaker form of the Doob- 
Meyer decomposition are however possible in the case of strong martingales. Namely, 
if M G M'g{T), then there exists a unique ^^"P^'^dictable increasing process [M](^) and 

there exists a unique .F^-predictable increasing process [M]^^) such that — [Mj^-* is an 
i-martingale for i = 1, 2. As noted in [3], for a strong martingale M, either [M]^^) or [M]^^) 
can serve as the process (M) above, but the question remains about whether the equality 
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cases, the answer to the latter is in fact affirmative. For example, if M G Al|(T) n A<^(T), 
or if M G A4| and (JF^) is a filtration generated by a standard two-parameter Wiener pro- 
cess (Wiener sheet), then [M](y = [M]^^) a.s. 

Although it is possible to develop the theory of stochastic integration in the plane 
with respect to general two-parameter martingales M € Al^(T) and define corresponding 
stochastic integrals J (pdM and J f ipdMdM (see [3] , [4] ) , as well as the so-called mixed area 
integrals (as in [4]) of the form J f hdjjidM and / f gdMdjj, (where /x^, ^ G T, is a continuous 
random function of bounded variation adapted to {Tz)-, and such that 1/^1(1") < C a.s. for 
some constant C < cxo, where |/i| denotes the total variation measure corresponding to 
the signed measure that n generates), but for the purposes of the present paper it will be 
sufficient to study such integrals in the special case when M is a standard two-parameter 
Wiener process (i.e. a standard two-paramctcr Wiener sheet). 

Recall that if is a random measure in , which assigns to each Borel set A a Gaus- 
sian random variable of mean zero and variance \{A), where A is the 2-dim Lebesgue mea- 
sure, and which assigns independent random variables to disjoint sets, then the stochastic 
process W = {Wz, z £ M^) defined by Wz := W{Rz), where Rz := (0,2;] is the rectan- 
gle whose lower left-hand corner is the origin and whose upper right-hand corner is z, is 
called a two-parameter Wiener process or a Wiener sheet. Equivalently, one could define 
a two-parameter Wiener sheet {Wz, z G M^) as a continuous Gaussian random field on 
with mean and the covariance function given by: 



Let {Wz,Tz,z G T} be a Wiener sheet. Let us introduce the following classes of 
integrands. Let {0^, z G T} be a process such that the following conditions hold: 

(a) (f) is a bimeasurable function of {u, z), 

(b) Jrj^Ecpldz < oo, 
and for each z G T, 

either (cq) (pz is ^z-measurable, 
or (ci) (l)z is J^l-mediSmahle, 
or (C2) (t)z is .?v^-measurable. 

Definition 1.4 For i = 0, 1,2, let Tii denote the space of (f) satisfying (a),(h) and (ci). 

Then one can show that for (f) G Hi, i = 0,1, 2, the stochastic integral (pzdWz can be 
constructed (as in [5]). Moreover, if one defines the process 



then the process (t)oW \s & strong martingale for ^ G Ho, a 1-martingale for cf) E Hi and 
a 2-martingale for (f) eH2- Moreover, define a process 



E{WzWz') = min(zi, 4) min(z2, 4)> Vz, z' G 



{<P oW)z= [ 4>cdW^ = I z)ct)^dW^ 



z G T 
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Then ^ = {^z,z G T) is a martingale with respect to {J^z)zeT if <?!',V' ^ 'Hoi a 1-martingale 
if ^Z" G Hi and a 2-martingale if ^, G ^2- In all cases continuous versions of the above 
defined processes can be chosen. 



Definition 1.5 Let H denote the space of functions ip{uj,z,z') = ipz,z'{''^) on J7 x T x T 
which satisfy the following conditions: 

(a) : is a measurable process and for all z, z' G T, ^z,z' is J^zvz' -i^^o,surable, and 

(b) : J J I{zX z')E{^l^,}dzdz' < oo. 

Then for arbitrary ip eH, the stochastic integrals 

JjRz-xRz JJr^xRz JJr^xRz 

are well-defined (as in [5]) for all z G T and X, Y^, Y'^ are respectively a martingale, an 
(adapted) 1-martingale and an (adapted) 2-martingale, and in all the cases the sample- 
continuous versions can be chosen. (See [5] for definitions of adapted 1- and 2-martingales.) 
Note also that the above double integrals are defined in such a way that only the values 
of the integrand on z X z' have an effect on each integral. 



Finally, the following proposition will be useful to us later on. 



Proposition 1.1 [6j Let [Xz,J^z',z G T = (0, Ti] x (0, r2]} be a strong martingale in 
Al|(T) satisfying either of the following conditions: (i) {J^z) is a filtration generated by a 
Brownian sheet; or (ii) X E A1^(T). Then 
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exp{Xz--{X), 



is a martingale iff E 
where T = (Ti, Tg) G M^. as before. 



exp{XT--(X)r} 



1, 



1.2 Fractional calculus and properties of fractional Brownian sheet 
Definition 1.6 Let ip{x) G Li(a, 6) (a,b eR) and let a > 0. The integrals 

(/.»W:=i^f .>«, (4) 

(MM:=j^/(^^*, .<a, (5) 

are called left-sided and right-sided Riemann-Liouville fractional integrals of order a. 

Fractional integrals (4) and (5) can, in fact, be defined for functions ip{x) G Li{a,b), 
existing almost everywhere. The following formula for fractional integration by parts is 
valid and will be useful (see [7]): For arbitrary ip{x) G Lp{a,b) and ip{x) G Lq{a,b), where 
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either p ^ + q ^ < 1 + a and (p,q) G N \ {(1, 1)}, or p = g = 1 but p ^ + q ^ < 1 + a, the 
following relation holds: 



j\{x){l2^i^){x)dx = f ^{x){I^_^){x)dx. (6) 

J a J a 

Moreover, fractional integration has the following semigroup property: 

/„VaV = ^atV /?-4^_^ = ^rV Va,/3>0, (7) 

where the above equation holds for every point in (a, 6) if G C[a, h] and for almost all 
points in (a, 6) if G Li{a,b). 

Definition 1.7 For functions f(x) on interval [a, 6] C M, i/te expressions (if they exist) 

where < a < 1, are called respectively the left-handed and right-handed fractional 
Riemann-Liouville derivatives of order a. Moreover, for a > 1, let [a] and {a} de- 
note, respectively, the integral part and the "fractional" part of a, < {a} < 1, so that 
a = [a] + {a}. Then the expressions (if they exist) 



i (±Y r fit) 

r(n -a)\dxj J a {x- ' 



:= ( - ) / t::— with n=[a] + 1, (10) 



(P« /)(a;) := J f [ f^^^ dt, with n=\a] + 1, (11) 

are i/te corresponding fractional derivatives of arbitrary order a > 1. 

For a < 0, we will also use the notation {I^_^(p){x) := (V^_^ip)(x) and {I^_ip)(x) := 
{V^^(p){x). Also define /°_,_ and /°_ to be the identity operators: Ia+(p = (f and I^_ip = (f. 

Definition 1.8 For a > 0, let I^_^{Lp) and I^_{Lq) be defined as the spaces of functions 
f(x) and g{x), respectively, of the form: 

f = I^+'P for some (f G Lp{a,b), 1 < p < oo, and (12) 

g = /"_V' for some V ^ Lq{a,b), 1 < q < oo. (13) 

Then fractional integration and differentiation are reciprocal operations in the following 

sense: For a > and arbitrary ip G Li{a, b), {'D^^I2^{p){x) = <p{x) for almost all x G [a, b], 
while {I2+'D^+f){x) = f{x) is satisfied for / G I"^{Li). Note, however, that there exist 
functions / ^ /^_|_(Li), whose fractional derivatives V^^f exist. 
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Definition 1.9 Let X be a finite interval. The function f{x) given on X is said to satisfy 
the Holder condition of order A G (0, 1] on X (or be Holder continuous of order X) if there 
exists a constant C > 0, such that 

\f{x)-f{y)\<C\x-y\^ (14) 

for all x,y G X. We denote by M.-'^{X), the space of functions satisfying (14). 
More generally, let X be a finite d-dim rectangle. The function f of x = (xi, . . . ,0;^) G X 
is said to satisfy the Holder condition of order (Ai, . . . , A^) G (0, 1]*^ on X (or be Holder 
continuous of order (Ai, . . . , A^) ) if there exists a constant C > 0, such that 

\f{x) - f{y)\ < C - yil^i + . . . + - yal^'^) (15) 

for all x,y E X. We denote by H'^i'''^<*(^), the space of functions satisfying (15). 
The following theorem will be useful to us later in the paper: 

Theorem 1.2 [7] Let f{x) = {x - a)-i^g{x), where g{x) G M.^{[a,b]), [a, 5] C M, A > a, 
-a < /X < 1. Then, f{x) G Ia+{Lp) if n + a < ^ for I < p < oo. 

Next let us focus our attention on the properties of a fractional Brownian sheet, which 
will serve as a model for the multiparameter observation noise in the nonlinear filtering 
problem discussed in Section 2. The interest in studying this type of random field stems 
from the fact that it has a number of remarkable properties which make it both mathemat- 
ically and practically interesting object, which is potentially useful in a large number of 
real-life applications. In fact, its one-parameter version, called fractional Brownian motion, 
has recently become an important modelling tool in geophysical and biophysical sciences, 
internet traffic modelling, financial applications and environmental sciences. The main 
properties of the fractional Brownian motion (fBm) are self-similarity, ability to model 
both short and long-memory effects (depending on the value of its Hurst parameter) and 
its non-semimartingale and non-Markovian structure. While the latter properties often 
make stochastic analysis of the dynamics driven by fBm very challenging, some established 
close connections with the standard Wiener process through fractional calculus techniques 
(some of which were mentioned earlier in this section) provide a number of mathematical 
tools to make it more tractable. 

Recall that a fractional Brownian motion (fBm) with Hurst parameter G (0, 1) is a 
continuous mean zero Gaussian process {B^ , t G M_|_), starting at almost surely, whose 
covariance structure is given by: 

^H{s,t) := E(Sf ) = i(|s|2^ + \t\^" -\t- s|2^), \fs,t G R+. (16) 

When H = ^, B^^ reduces to a standard Wiener process (or standard Brownian motion). 
For > i, the increments of the fBm are positively correlated and the fBm exhibits 
long range dependence (long-memory property): Yll!Li \{Bi — Bq){B^j^^ — B^)\ = oo. 
When < i, the increments of the fBm become negatively correlated, resulting in its short 
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memory. For every H G (0, 1), the £Bm is a self-similar process with self-similarity 

index iJ, since {B^)t>o = {c^Bf^)t>o for all constant c > (where = denotes equality of 
two processes in distribution). Also, clearly, for any H ^ ^, is not a semimartingale, 
implying that the standard techniques of stochastic calculus and stochastic integration 
are not directly applicable in the fBm case. Moreover, sample paths of B^ are nowhere 
differentiable (with probability one) but the trajectories are Holder continuous of any order 
strictly less than H. Finally let us note that the fBm enjoys a number of fractional integral 
relations with respect to a standard Wiener process. For example, B^ = Knit, s)dWs 
for some standard Wiener process W, where 

Knit, s) = (( J)"""" {t - sf-'^ -{H- \)sh-^ u"-l {u - sf-'^du^ , (17) 

/ 2HV( — H) 

where ch = \ p,„ It is useful to note that Knit, s) can also be represented 

by: 

KH{t,s) = c^s^^(lf"^u^-h[o,i](n))(s), (18) 

H-- 

where = ch^{H ^) and ^ is the right-sided Riemann-Liouville fractional integral 
of order H — 1/2, introduced earlier in this section. The above Wiener process W can be 
reconstructed from B^ via Wt = Jq K^^{t, s)dB^ , where the kernel K^^^ is given by 



K^Ht, s) = 4 ( Q ''~ht-s)'^-"-{H-l)s'^-« {u - s^^-^dn^ , 



(19) 



where c'rr = ^, . / ^^.^.a*-^^^^,-*^ i, . The latter kernel can also be written in the fol- 
lowing form: 

KH\t,s) = ^.^^(/i-''u^-h[o,](n))(s). (20) 

Note that the above processes B^ and W generate the same natural filtrations. 

The two-parameter fractional Brownian sheet (fBs) with Hurst indices {a, f3) € (0, 1)^, 
represents the two-parameter analogue of fBm, and can be defined as a continuous centered 
Gaussian random field B°''/^ = {Bz'^,z G K^), whose covariance structure is given by: 
E(i3"''^i?"'^) = 'jaizi, z[)-jj3{z2, , with 7q,7/3 defined as in (16) and \fz = {zi,Z2), z' = 
{z'l, Z2) € M^. Naturally the fBs inherits all the remarkable properties of an fBm, while 
allowing one to introduce some new effects (like having long memory in one parameter, and 
short memory in the other, for example). In the case when both Hurst indices are greater 
than 1/2, we will call the fBs persistent, as it displays long memory in both parameters. 
Note also that the fBs could also be equivalently defined through its integral representation 
with respect to a standard Wiener sheet. Namely, 

B^'^:= JJ^ KUzi,Ci)Kpiz2,C2)dW^<:,X2)^ z e^l, (21) 
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where Rz denotes the rectangle (0, z] = (0, Zi] x (0, Z2] (as before) and where (W^)^gR2^ is a 
standard Wiener sheet. On the other hand, if we let K~^p{z;Q = K~^{ziXi)Kp^{z2;C2), 
where K~^{s,t) is the kernel defined in (19) (or (20)), then the following representation 
is valid: Wz = K^^g^z; ()dB'^'^ . Clearly, above definition and properties extend to pa- 
rameter spaces of dimension higher than two, but for brevity we will restrict our attention 
to the two-parameter fBs case throughout the paper. 



2 Nonlinear filtering of random fields with persistent frac- 
tional Brownian sheet observation noise 

While in the one-parameter case the topic of optimal nonlinear filtering with fractional 
Gaussian observation noise has been studied quite extensively in a variety of contexts (see 
e.g., [8], [9], [10], [11], [12], [13]), there is currently no mathematical literature devoted 
to similar questions in the context of "spatial" filtering of multiparameter random fields. 
Thus, we expect that results presented in this paper will be of interest to both theoretical 
and applied scientists, especially in view of an increasing use of imaging technology in a 
number of fields (ranging from biomedical applications to surveillance), where the spatial 
structure of the underlying "signal" of interest is important and denoising and filtering of 
"noisy" images and videostreams are clearly needed. 



2.1 Observation model with persistent fBs noise. "Fractional-spatial" 
B ayes' formula. 

Consider the following observation model: 

Yz= f g{X^)dC + B^'f^, zeT=[0, Ti] x [0, T2] C M^, (22) 

where Rz = (0, z] = (0, zi] x (0, Z2], the signal of interest X = {Xz, z &T) and the obser- 
vation random field Y = {Yz,z G T) are measurable (.?vj:)-adapted random fields defined 
on a complete filtered probability space {^,J-, {J-'z),P), where the filtration (J^z) satisfies 
conditions {Fl)-{F4) given in Section 1.1, and B'^'^ = {Bz'^ , z G T) is a fractional Brow- 
nian sheet on (^i,^, {!Fz),P) with Hurst parameters a,(5 ^ (^, 1) and B"^'^ is assumed to 
be independent of the signal process X. Throughout the paper let us assume that the 
following conditions are satisfied: 



(Ai) Function 51 : M ^ M is Holder-continuous of order A on any finite interval in 
where A > 2max(a,/3) — 1; 
and 

(A2) The following integr ability condition is satisfied: 



S2 



E 



'2?0+V(^))(Cl,C2) 



<^ClC^C2 < 00, 



(23) 
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where ^ , Vq_^_ ^ are the fractional Riemann-Liouville derivatives defined in Defini- 
tion 1.7, 

gtiX){u;) := zr'"zr'^giX,{u;)), Vz = (^1,^2) G T, Vu; G Q, (24) 

and "(2>" denotes the tensor product of operators. Namely, given a function / : T — M 
and a pair of linear operators £±,£2, defined on appropriate functions of the form /i : 
[0,ri] R and /2 : [0,T2] ^ M, respectively, let 

A®/:2/j(2i,^2) := A(>C2/(-,^2))(2i), v(zi,z2) € [o,ri] X [o,r2]. 

Lemma 2.1 Fix arbitrary a,l3E (^, l). Let h : T ^ M. be a Holder continuous function 
of order (Ai, A2), where Ai > ct—^ and A2 > Then there exists a function (5/i : T — > R 

such that 5h G L2(T) anc/ 

/ K-Uzi,Z2;0h{0dC= f 5„(CK, V(zi,Z2)GT. (25) 

^[0,Zl]x[0,22] V[0,2i]x[0,Z2] 

Moreover, if we let h*{zi,Z2) ■= z^ Z2 h{zi,Z2), then Sh can be taken as follows: 
Sh{zi,Z2) = ^zr'^ztHv^+~'^V^+~'h*)izi,Z2), y{zi,Z2) e [0,Ti] X [0,T2]. (26) 

Proof: Since linear combinations of tensor products of functions of single variable are 
dense in the space of functions of two variables, it suffices to consider the case of functions 
/i of the form h{zi,Z2) = /ii(^i)/i2(^2), where hi G H^i([0,ri]), /12 G M^^{[0,T2]). By (20), 

/ K-^^izi,Z2;0h{0dC = ^ rCr"«"Jn-h[o,,,](u))(Ci)/ii(Ci)dCi 

J[0,zi]x[0,Z2] ^a^/i Jo 

r,_l , 

X / Cr\K~-^^~~'M0M(^))(C2)h2{C2)dC2 

Jo 



cr ^v;--'u-^--hl{u)){a)dCl 

Sh{Ci,C2)dCidC2, 



[0,Z1]X[0,Z2] 



^\l^'{V^0^^U^-^h2{u)){C2)dC2 



where 5h is defined by (26) and where wc used the fractional differentiation by parts 
formula (see corollaries from (6) in [7]), together with Theorem 1.2. Note also that if 

hi G ll^i([0,ri]), where Ai > a - i, then, by Theorem 1.2, hi G {L2{[Q,Ti])), which 

implies that /ii(s)ds G /q^ ^ (L2([0, Ti])). If we let 



6hAs) := \s''-'2{v;^K^2-"hi{v)){s), ys G [o,ri], 



(27) 
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then one can easily check that 

K^{t,s)Sh,{s)ds= f hi{s)ds, ViG[0,Ti]. 
^0 

Since the integral operator /Cq associated with the kernel Ka, i.e. 

[/Ca/](t)= f K^{t,s)f{s)ds, /GL2([0,Ti]), 
Jo 

is an isomorphism from L2([0, Ti]) onto 1^^ ^ (L2([0, Ti])), then / G L2([0,ri]) if and only 

if ICaf € /o+^(L2([0,ri])). Therefore, 4^ G L2([0,ri]). Similarly one shows that if 
/i2 G H^2Qo,r2]), where A2 > /3 - 5, then 6^2 ^ ^2([0,r2]), where 

Sh^s) := Vs e [0,Ti]. 

Therefore, 5hi®h2 = (^hi ® (J^z S L2([0,Ti] x [0,T2]), and the required result follows. □ 

Corollary 2.2 Fix Aq G ( ™'^^*'^^^ — Suppose the signal X = G T) has almost 

surely Holder-continuous sample paths of order (Ao,Ao) and g satisfies condition (Ai). 
Then for almost all to E ^ one can define function {Sz{X),z G T) by 

S,{X)iu) := -^z^''zt^V^-'^^V^,;'^g:iX)iu,)){z), z = (^1,^2) G T, (28) 

with g*{X) defined by (24). Then, assuming also that (A2) holds, d{X) = {dz{X),z G T) 
has the following properties: 

(i) S.{X){ui) G L2(T) for almost all uj e and E J^{5z{X))'^dz < 00; 

(a) For every rectangle Rz = [0, C T, 

I K-J(z;05(^cK= / (^cWC a.s. (29) 

J Rz ^ Rz 

From now on suppose that the assumptions of Corollary 2.2 are satisfied. Let us 
introduce processes 

Wl:=[ K-],{z-OdY^ ^nd := f K-],{z- Qdsf , z e T . (30) 
Then it is easy to see that = Jj^J^{X)d( + VFf . Next let us define a process 

v = {Vz,zeT)hr. 

F, = exp|-/' 5^{X)dW^-]-j {5^{X)fdc\, zeT. (31) 

k J Rz *f Rz ) 

Note that 

y, = exp|-/" 5^{X)dWl + ]-l {5^{X)fdcV zeT, 
K J Rz J Rz ) 



12 



thus, 

y, = exp|-^ SdX)d{J^ K^^{C-C')dY^,) + \j^ {5^{X)fdC^, z G T. (32) 
Lemma 2.3 Let V = {Vz,z € T) be defined by (32) (or, equivalently, by (SI)). Then 

K(%i,T,)) = 1. 

Proof: Since B"'^ and X arc independent, then and X are independent, which 
imphes that one can define a standard Wiener sheet on a complete probabihty space 
{^2,^2, P2), define X on a complete probability space (J7i,^i,Pi) and then consider the 
processes on a product probability space (Jli x 0,2,^1 x ^2, Pi x -^2)5 with W^{lo) = 
W^{uJ2) and X{uj) = X{uJi) for all lo = {uji,uj2) G Jli x ^2- Let Z{uj) = Z{uji,uj2) = 
J^6^iX{uJi))dWf{uj2). Then, upon taking into account Corollary 2.2, it follows that for 
almost all (fixed) uoi, Z{ijJi, •) is a Gaussian random variable with mean and variance 
Jrj, [di^{X{uJi))]'^d(^, and, thus, for almost all wi G fli, 

Ep.(%,T.)(^i,-))=Ep, exp|-Z(a;i,-)-^[^c(^(^i))]'<} =1' 
which implies that 

E(%i,T2))= / F(Ti,T2)(t^i,'^2)(Pi xP2)(da;i,dw2)= / lPi(da;i) = l. □ 

Theorem 2.4 Consider observation model (22), where the signal X = {Xz,z G T) has 
almost surely Holder- continuous sample paths of order (Aq, Aq), where — 2. < < 
^, and g satisfies conditions (Pl\)-(K2). Let P be a new probability measure on {0,,!F) 
given by: 

dp 

— = V^T„T,) a.s.{P) (33) 

Then P is equivalent to P and, under P, (22) holds a.s., Y is a standard fBs with Hurst 
indices {a, (3), X has the same law as under P, and processes X and Y are independent 
under P. Moreover, the following "spatial-fractional" version of the Bayes' formula holds: 
For any F G Cb{R), 

E[F{X,)V,-] ATj] E\F{X,)V-^\J^Y] 

where E denotes the mathematical expectation under P, J-'J denotes the filtration generated 
by the observation process in the rectangle Rz = [0, zi] x [0,^2] C T, i.e. 

■.= a{Y^ :0^C^2), -zeT, 
and V = {Vz,z G T) is defined by (32) (in terms of (28)). 
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Proof: The first part of the theorem follows at once from the multiparameter Girsanov- 
type theorem for the standard Wiener sheet (see e.g. Theorem 1 in [15], p. 89), Proposi- 
tion 1.1 and Lemma 2.3. To prove Bayes' formula, arguments similar to those constructed 
in [2] for the two-parameter Wiener sheet observation noise can be used. The proof, to a 
large extent, follows the lines of standard arguments used in the one-parameter (fractional 
noise) case, thus, we will omit the details here. □ 

2.2 Evolution equation for the optimal nonlinear filter along an arbi- 
trary increasing 1-dim curve 

Here we present a stochastic evolution equation satisfied by the unnormalized optimal 
filter when its dynamics is tracked along an arbitrary monotone non-decreasing 1-dim 
continuous curve A connecting the origin to the point T = {T\,T2). By a monotone 
non-decreasing path we mean that A is nondecreasing (in the sense of partial ordering 
in the plane) in both zi and Z2 directions. For each z G T = [0, Ti] x [0, T2], let be 
the "smallest" point on A which is larger than or equal to z with respect to the partial 
ordering ^. The path A divides domain T into two regions; the region below A, which is 
denoted by and the region above A, denoted by 

Namely, = {( g T : COCa = Ca} and = {C G T : Ca ©C = Ca}, where we use the 
notation a 6 := (ai, 62) for arbitrary a = (ai, 02), & = (&i, ^^2) G 1^+ (as in Section 1.1). 

Definition 2.1 Let [TziZ G T) 6e a filtration satisfying conditions (F1)-(F4) of Sec- 
tion 1.1. Suppose A is a monotone nondecreasing continuous 1-dim curve connecting the 
origin to point T = (Ti,T2) G M^. Then 

i) A process (j) = {(j)z, z gT) is called A-adapted if (f>z is J- z ^-measurable for all z G T. 
a) A process X = {Xz,z G T) is called a A-martingale if X is A-adapted and 

E [X{z, z'\ I =0 for all O^z^z' ^ T. 

Definition 2.2 Let Ha be the space of processes cj) = {4'z,z G T) satisfying the following 
conditions: 

(a) (j) is a bimeasurable function of {io,z); 

(b) JjE(j)ldz < 00; 
(ca) (p 'i'S A-adapted. 

For G Ha, define processes <pf^ = {4>f^,z G T) G Tij, i = 1,2 (see Definition 1.4 for 
definitions of 'Hi,TC2), by: 



' 0, otherwise; 1 0, otherwise. 



Then <f>z = + (t>2z almost all z G T and one can construct stochastic integral 
Jrj. (f>zdWz = {(j) o W)^ for 4> G Wa and show the following properties for the resulting 
integral (see [14] for details): 
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Proposition 2.5 Let A be a monotone nondecreasing 1-dim continuous curve connecting 
the origin to the final point T . Let (f) G TYa md define the stochastic integral of 4> with 
respect to a standard Wiener sheet {Wz,J^zjZ G T) by 

((/.oT^)^ = (</,f oT^), + (0^oT^)„ zeT, (35) 

where the two stochastic integrals on the right-hand side of (35) were discussed earlier in 
Section 1.1. Then the integral has the following properties: 
i) ((f> o W)'^ is a A-martingale; 

a) {(f) o W)^ is a one-parameter martingale on the path A; 

Hi) If A and A' are two monotone nondecreasing paths connecting the origin to T and both 
passing through a point zq G T, and^ is both A and A' -adapted, then {(poW)^^ = (^oVF)^ . 

Suppose our signal Xz is a two-parameter semimartingale in the plane of the form: 

Xz = Xo+ j (l)^dW^+ j 0^d(:+ ! (36) 

+ f fu'dCdWc + I gcC'dWcdC, z e T, (37) 

J Rz ^ Rz " Rz ^ Rz 

where, as usual, Rz = [0, zi] x [0, Z2] and 4> eHo and ip,f,g G Ti, where spaces T~Cq,T~C are 
defined as in Definition 1.4 and Definition 1.5. 

Then, by [14], for an arbitrary monotone nondecreasing continuous 1-dim curve A, 
connecting the origin to T, there exist rf(^ = r){A, Q and z/^ = iy{A, Q such that 77 G Ha 
and 

Xz = Xo+ [ 7i{A,0dW^+ [ i^{A,Od(:, z€A. (38) 
•JRz -IRz 

If 6 is A-adapted, then v can be chosen A-adapted. As such, X is clearly a sample- 
continuous semimartingale on A. 

In the rest of Section 2.2 we will therefore assume that the signal process X is of the 
form (38) , where the standard Wiener sheet W is independent of the observation random 
field Y . Let us consider the nonlinear filtering model (22) along with conditions (A^), (A2) 
and recall the general framework of Section 2.1. 

Theorem 2.6 Let A he an arbitrary monotone nondecreasing continuous 1-dim curve 
connecting the origin to the final point T G M^. Let us assume that the observation 
model (22) holds, along with conditions (A-^), (A2), and suppose that the signal X is a two- 
parameter semimartingale in the plane, which is written in the form (38) , where 77 G T^a 
and V is A-adapted, and whose trajectories are Holder- continuous of order (Aq, Aq), where 
Ao > ^^^"'^^^ 2 _ Pqj. p g C^(M), consider the unnormalized optimal filter 

az{F):=E[F{Xz)V-^\Tj], zeT, (39) 

introduced in Theorem 2.4- Then the following stochastic evolution equation, governing 
the dynamics of the unnormalized optimal filter along the monotone increasing path A, is 
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satisfied: 

a,{F) = ao{F) + l^ a^^{„F' + ^r,^F")dC + i^-^(C; CO^^^C') > ^ A, 

(40) 

where 

:= E [F {X^J S^^{X)V^l^ \ , (41) 
and {Sz{X),z € T) is defined in (28). 

Proof: Let us reparamctcrizc A by {z{t); < t < 1} so that the process {Xz, z G A} can 
be rewritten as {X^^^^jO < t < 1}. By Proposition 2.5, X is a continuous one-parameter 
semimartingale on A, thus, by Ito's formula (for one-parameter case), for all F G C^(M), 

F(X,(,)) = F(X,(o)) +^V'(X,(,))dX,(,) + i^*F"(X,(,))d(X,X),(,), t G [0,1], 

where (X, X)z(^t) = /r 'V^'^C- Note that one can re-express F along A free of the earlier 
parametrization as follows: 

F{Xz) = F{Xo) + [ F'{X^JdX^ + ^[ F"{X^^)rildC. z e A. 

Similarly, since 

=exp|^^.(.)(^)dW-,t,) + ^_^*[<5.(.)(X)pciz(s)|, tG [0,1], 

then 

= ^ + jyz~ilM)iX)dKs) + jy^'il)^^^^^^ *e [0,1], 

where the latter equation can also be rewritten free of parametrization as 

V-^ = l+[ V^l^S^^{X)dWl, zeA. 

Moreover, 

Ki)nX.^t)) = F{X,^o)) + + ^r?,%)F"(X,(,))) dz{s) 

+ jy.~i^)F\Xz^,))r,,^s)dW,^s) + jyiX,^s))V^^^^ t G [0,1]. 

Then, upon taking conditional expectations of both sides of the above one-parameter 
equation with respect to .F^^^ = JF^-j under P, one arrives at the following equation 
along the path A: 

az(t){F) = ^z{0){F) + f^az^s){^F' + \rfF'')dz{s) + j\ t € [0,1], 
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where 

The latter evolution along the 1-dimensional path A can be expressed free of parametriza- 
tion as follows: 

a,{F) = ao{F) + [ a^JuF' + \ifF")dC + / a^^{F5)dWl , z G A, 
or, equivalently, 

a,{F) = ao{F) + J a^^{vF' + \rfF")dC + j a^^{F5)d{^ J i^-J(C; CO^^^') > ^ € A, 

(42) 

where the equations hold almost surely under P and P. □ 

Note 2.1 Let us observe that, in contrast to the case of filtering in the presence of a 
martingale observation noise, the above stochastic evolution equation (40) cannot be in- 
terpreted as a measure-valued SPDE in view of the special meaning assigned to a{ ) in (41). 
The latter is necessary because 5z{X) is not a function of Xz but rather is a function of 
the entire "history" (-^^,0 -<C, < z). 

2.3 Analogue of Duncan-Mortensen-Zakai equation for the optimal filter 
in the case of 2-parameter dynamics with fractional Brownian sheet 
noise 

The stochastic evolution equation (40) developed in Section 2.2, governing the dynamics 
of the unnormalized optimal filter, is two-dimensional in form, but clearly one-dimensional 
in spirit. Our objective in this section is to develop a "fractional-spatial" analogue of the 
Duncan-Mortensen-Zakai equation for the unnormalized optimal filter which is inherently 
two-dimensional . 

Let a : M ^ R and b : M ^ M be measurable functions satisfying the following Lipshitz 
and growth conditions: there exists a finite constant C > such that for all x, y G M, 

\a{x)-a{y)\ + \h{x)-h{y)\<C\x-y\ 

and 

\a{x)\ + \h{x)\<C{l + \x\). 

Then there exists a unique strong solution to the following multiparameter SDE (see 
e.g. [15]): 

Xz = Xo+ I a(X^)dC+ j h{X^)dW^, zGT, 

J Rz J Rz 

where W denotes a standard Wiener sheet. Moreover, the solution has Holder-continuous 
sample path of order (Ai, A2) for all Ai, A2 G (O, \). 
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Theorem 2.7 In the framework of Section 2.1, assume that the observation model (22) 

holds, and that above Lipshitz and growth conditions on o(-) and b(-) are satisfied and 
conditions (A^^), (A2) are valid. Suppose that the signal X is the unique strong solution of 
the following SDE: 

X^ = Xq+ I a(X^)dC+ ! h{X(^)dW^, zgT, (43) 

J Rz J Rz 

where W is a standard Wiener sheet independent of the observation Y. Let a'z{F) := 
M[F(Xz)V~^ \ i.e. CTz(F) is the unnormalized conditional expectation corresponding 

to the optimal filter. Then for all F G C^(R), evolution of the unnormalized optimal filter 
has the following structure: 

oz {F) = ao (F) + / ac (aF' + i b^F") dC + / {FS) dWY 
JRz ^ JRz 

+ // (Jc^^c,,{F;b®b)dWldWl, 

JJRzXRz 

JjRzXRz ^ 

+ [[ [cTcciF'- S®a) + la^,c'{F"; 6 b^)] dW^dC 
JJRzXRz ^ 



+ [[ liCxC') \ac,c'{F"; am) + l-a^,c{F"'; b2®a+a®b2) + ^ac,c'(F(-); b^^b^) 
JjRzXRz L ^ 4 ^ 



dCdC, 



(44) 

where denotes the tensor product of functions, az{F5) := E.[F{Xz)5z{X)V,, 
(^z,z'{F;5 ® 6) := E[F{Xzvz')Sz{X)Sz'{X)V;;;^],, l^Ivz']' for arbitrary functions fi : 

R ^ R, /2 : ^ M, we put az,z'{fi; f2) := E[/i(X,v.')/2(^., ^.OKvl' 1 -^w.'] M all 
z,z' G T. (In (44), Wj = /^^ K~'^p{z;C,)dYQ and 5 is given by (28), as before.) 

Note 2.2 In Theorem 2.7, we could write cFz{F) = az,z{F;l), where 1 denotes function 
on which is identically equal to one. 

Proof of Theorem 2.7: First note that, under P,Y is a fraetional Brownian sheet with 
Hurst indices (a,/3), while the corresponding field , given by = Jj^ K~^p(z; C)dYi^, 
is a standard Wiener sheet and the two random fields generate the same natural filtration, 
thus, the observation sigma-field {TY)q^z^t has properties (F1)-(F4) of Section 1.1. 
Similarly, {J'^)zeT ^-iid {J-^'''^)zeT have properties (F1)-(F4) under reference probability 
measure P. Note also that the paths oi X = {Xz,z £ T) are almost surely Holder- 
continuous of arbitrary order (Ai,A2), where Ai,A2 < ^, thus S is well-defined and the 
conclusions of Corollary 2.2 hold. Next, by a version of the Ito's formula for multiparameter 
semimartingales (see [5]), we obtain that for arbitrary F e 



F{Xz) = F{Xo) + / F'(X^) [aiX^)dC + biX^)dWc] +\ j F" {X^)h\X^)dC 
J Rz J Rz 
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+ 



+ 



+ 



RzxRz 



[F"(X^vf')K^c)a(Xc') + ^F'"{X^^^,)biX^)b\X^,)]dCdW^, 



■ / / [^"(^CvC')K^c)'i(^C') + lF"'{X^^^,)b{X^)b\X^>)]dW^dC' 

+ [[ i{c^-C')\f"{x^vcO<x^Mx^') + If'"{x^^<:0{<x^^^^ 

JJR^xR^ L ^ 

dCdC'. 

Similarly, under P, one shows that 

= 1 + / V^'^5^{X)dWl V^^\, 5(^{X)5(^,{X)dWldWl. 



+\F^^^\x^^^,)b\X^)b\X^,] 



a.s. 



Then the multiparameter version of the stochastic integration-by-parts formula (together 
with independence of W and under P) yields a corresponding equation for the product 
F{Xz)'V^^ . Upon taking conditional expectation of both sides of the latter equation for 
F{Xz)'V^'^ with respect to (where note that = TY ) and using Lemma 2.8, which 
is proved below, one arrives at the following equation: 

E(F(X,)y-i|J-r)=E(F(Xo)|^o'")+^ E(^[a(Xc)F'(Xc)+^b2(Xc)F"(Xc)]F^-^|^j)dC 

+ / t{F{X^)6c^{X)V^^\Tl)dWl 

J Rz 

+ If ^F{X^y^,)d^{X)dc,{X)V^-\, 

+ E([a(^c)^'(^CvC') + ^t''(^c)^"(^CvC')]^C'^CvC' 
■ / / K([a{X^,)F'iX^y^>) + ib2(X^,)^"(^cvC')]<5cVc' 



+ 



+ 11 /(CAC')E 
'RzXR^ 



F"(Xcvc')a(^c)a(^C') + ]F(^''\X^^^,)b\X^)b\X^,) 



+-F'"(X^vc'){^''(^C')a(Xc) + a{X^,)b\XO} 
thus, the required conclusion follows. □ 



^C^C'^dCdC a.s., 
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Lemma 2.8 Let W and be independent standard Wiener sheets on a probability space 
{n,J^T,P) and J^V'^"" := a{W^,W^ : ^ C ^ z,0 ^ C ^ z), z e T. Also let {J^f) and 
) denote the natural filtrations generated by W and , respectively. Consider a 
process M (which is [j^Y'^ ^-measurable), given by 

M, := [ (j)^dWl + [[ i^^^^'dW^dW^, 
Jrz JJrzxRz 



where (p € Ho and ijj ^H, with 7^o and H being defined with respect to filtration {T-, 
Then 

i) For any process ip Eli, 



E 



E 



RzxRz 



a.s. P, 



J=T 1 = a.s. P, 



\JRz->^Rz 

E ( / VCC'^^C^^C' 

\JRzXRz 



E 



RzxRz 



a.s. P, 
= a.s. P, 
= a.s. P. 



ii) The following equation holds almost surely with respect to P: 

E (m, I J^-f ") = f E (</>J ") dWl + [ E (Vcc I ^c'^c") dW^dW^. 
*t Rz J Rz X Rz 

Also, 



E 



Y 



RzxRz 



^p<:,cdCdw^', 



^ a.s. P, 



E 



(^1^ ^iX'dW^dC' 



'RzXRz 



= / E{i^c,C'\^^c')dCdW, 

/ JRzXRz 

HA,C'\^0/c)dWldC' a.s. P. 



RzXRz 



Proof: Let us start by showing that the first equahty in (i) holds, i.e. that 



E 



jRzXRz 



^<dW(^dW^'\J^, 







almost surely under P. By independence of W and , we may assume that is a 
standard Wiener sheet on a filtered complete probability space {i}^ ,{J^^)zeT,P^), 
whereas is a standard Wiener sheet on another filtered complete probability space 
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{DX ,!F^^ ,{Tf^^)z^T,P^), where {!F^) and (T^^) are, respectively, natural filtrations 
generated by processes W and (in and 17^, respectively), and {VL,Tt,P) = 

{Q^ xQ^ xj'^'' xP^), i.e. the product probability space. Then W and are 
defined on (O, Jr,^) by W^iu) = W,{uJi) and (uj) = Wj {uj2) for all oj = (wi.tua) G O. 
Then, clearly, JP-f = J^f x {0,J7^} and = {0,^7^} x . 

Next let us fix some n G N and consider a partition of rectangle Rt = (0, T] (where 
T = (Ti,T2)) into rectangles A^^ := (z(,,j), , where z^ij) = (2-"zTi, 2-'^ir2). 

Let (S be the class of processes ^p of the form: 

2"-l 

v^(c,c')= E «.iHiA„,(c)iA,.,(c'), (45) 

i,j,k/=0 

where aijke is ^l^^'^^vz^t «x -measurable. By definition of the double integral, 

„ 2"-l 

/ V(C,C')c;W^C'^^C' := I] aijkel{i<k}hi<j}W{R,nAij)W{R,nAk,i). 

•I Rz X iiz • ■ 7, «_n 



i,j,k,e=0 



Then, 



if / ^(CCO^^w^c^^c'l-^f'l 

oiiouhi<k-)^^{^<3}^{Pz n Ai,j)T^(i?^ n Afe,^) I {0, O^} X J', 



= E ^ 

Note that VQ G JF]^'', we have 

/ oiijki{uii,u}2)l{i<k}n{e<j}W{Rz n Aij)(iOi)W{Rz n Afe/)(u;i)(iP(c<;i,a;2) 

ajjfc^(u;i,u;2)l{i<fc}n{£<i}^(^z n Aij)(wi)M/(i2^ n Ak,e){oJi)dP^ (coi) dP^ {0J2) 

= 0, since for fixed uj2 G ^^2 and for all i < k and £ < j, random variables ajjfe^(-, 0^2), 
W{Rz n Afe,^) and VF(-R2 fl Ajj) are mutually independent. Thus, 



E 



a.s.(P), 



aijkil{i<k}n{e<j}W{Rz n Aij)W{Rz n A^,^) | {0,fi^} x j^f 
implying that for any simple process E S, 

fif / i^iCOdW^dW^'lJ'^''] = a.s. (P). 

\JRzXRz J 

Since 5 is dense in the required equality follows for arbitrary ^ G H by taking ap- 
propriate limits. Similar arguments show that the remaining equalities in (i) are also 
valid. 
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To prove (ii), let us show that VV' G 

^{ I ^{C,C,')dWldWl,\J^r\ = [ Eii^iCaiJ^^ljdWldW^ a.s. (46) 
First, consider ^ G iS of the form (45). Then 



E 



'RzXR^ 
2"-l 

= ^<^iM\^Y'']hi<k}n{i<j}W^{Rz^^i,j)W^{Rz^^M), (47) 



where note that 

tY 



^{i<k}n{t<3}'^ {Rz n \,j)W {R^ n Akt) = unless Z(^k,j) = V z^k,£}) z. 

Since a^jk^ is ^1^'^ -measurable and ^^(fej) z, then the conditional expectation in the 
right-hand side of (47) satisfies equation 

^aijkelJ^^''] = E[aijke\^^^^J a.s., 

by independence of W and and since Wiener sheets generate independently scattered 
measures. Thus, 



e( [ ^{C,C')dW^dW}\J-r) 

KJR^xR, J 



2"-l 



= ^.^im\^7^^u,i^\i<k}v^{i<,}W'' {R,^ Ai^^^^^ a.s. (48) 

On the other hand, 

/ E[^(C,Ol^o^J'']^^c'^^c 

JRzxRz 

2" — 1 

= / E[a,,fe,|J^^^j:]lA,,(C)lA,.,(C')^W^c'^^C 

2"-l 

E lo<fe}n{€<,}E[a^,-fc^|^^3lA,,(C)lA,,,(C0^^W-c'^iW^C, (49) 

'^^^ i,j,k,e=o 

where the last equality holds by definition of the double integral. Note that 

l{i<fe}n{^<j}lAi,,(C)lAfc,^(C') / implies that C V C' G A^^^^ , 
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which, in turn, implies that E[ayfe^| J^^^,] (in the right-hand side of (49)) equals almost 

surely to E [ajjfc^lJ^j^^^ J , since aij^e is -^i^^ -measurable. Thus, from (49) by definition 
of the double integral, 

/ ^^p{c,c')\:F^l]dwldw^ 

JRzXRz 

= ^<^iM\^YZ)]hi<k}n{t<j}^^{Rz^^i,3W''iRz^^^^^ a.s. (50) 

From (48) and (50), it follows that (46) holds for all tfj G S. Since S is dense in Ti., it follows 
that (46) holds for all t/' € ?Y by taking appropriate limits. Similarly one establishes that 
V</. G Ho, 

E( f ^i;dWl\Tf''') = / E{^^\J^^'^)dWl a.s., 
</ Rz ^ Rz 

thus, the first statement in (ii) is proved. The remaining two statements in (ii) can be 
established by analogous arguments. □ 

3 Conclusions 

In this paper the problem of spatial nonlinear filtering of a multiparameter semimartin- 
gale random field, with estimation based on an observation random field perturbed by 
a long-memory fractional noise, has been considered. Two types of stochastic evolution 
equations, governing the dynamics of the unnormalized optimal filter in the 2-dimensional 
plane, has been derived. One equation follows the dynamics of the optimal filter along an 
arbitrary non-decreasing (in the sense of partial ordering) one-dimensional curve, while 
the other describes behavior of the optimal filter in terms of "truly" 2-dimensional dynam- 
ics. In view of long-memory in the observation noise, neither equation can be viewed as 
measure-valued SPDE and their interpretation is not trivial. However natural questions 
regarding uniqueness and robustness of the solutions to the evolution equations, as well 
as construction of suboptimal filters, can be addressed and the authors plan to do so in 
the forthcoming work. 

Despite numerous important practical applications of spatial nonlinear filtering (in 
connection with "denoising" and filtering of images and video-streams in physical, bio- 
logical and atmospheric sciences, for example), there currently appears to be very little 
mathematical literature on the subject. In particular, the results presented in this paper 
represent the first mathematical results pertaining to spatial nonlinear filtering of random 
fields in the presence of long-memory (fractional) spatial observation noise. 
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